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ABSTRACT 

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluc- 
tuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At 
small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become in- 
creasingly obscured by extragalactic foregounds and secondary CMB signals. We present results 
from a nine-parameter model describing these secondary effects, including the thermal and kinematic 
Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared 
Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extra- 
galactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of 
the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT 
data, fitting this model to the multi-frequency spectra in the multipole range 500 < £ < 10000. We 
extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 
220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model 
provides an excellent fit to both datasets simultaneously, with x 2 /dof= 675/697 for ACT, and 96/107 
for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT 
in bandpowers, marginalizing over the secondary parameters. This provides a simplified 'CMB-only' 
likelihood in the range 500 < I < 3500 for use in cosmological parameter estimation. 
Subject headings: cosmology: cosmic microwave background, cosmology: observations 



1. INTRODUCTION 

Measurements of the Cosmic Microwave Background 
(CMB) have played a central role in constraining cos- 
mological models. Anisotropies measured over the whole 
sky by WMAP have provided evidence for a flat uni- 
verse described by just six cosmological parameters. The 
measurement of the Sachs- Wolfe plateau in the power 
spectrum, and three acoustic peaks, have led to con- 
strai nts on ACDM parameters to percent-l evel accu- 
racy (|Komatsu et al.ll2011l:lLarson et al.ll201lD . The Silk 
damping tail of the power spectrum provides a wealth 
of additional information about the physics of the early 
universe, encoded in its shape, and in the p ositions an d 
heights of the higher-order acoustic peaks (jSilkl 119681 ) . 



Extracting information from these angular scales is com- 
plicated by the presence of additional power from ex- 
tragalactic point sources, emission from the Galaxy, and 
secondary anisotropies due to the thermal and kinemati c 
Sunyaev-Zel'dovich effects (jSunvaev fc Zerdovichlll970l) . 

The Atacama Cosmology Telescope (ACT) mapped 
the mm-wave sky with arcminute resolution from 2007 
to 2010 in two distinct areas. About 600 square de- 
grees were used to compute the angular power spec- 
trum. Power spectra and cosmological results us- 
ing the 1-year d ata, from the 2008 observing season , 
were presented inlFowler et al.l (|2010D ; iDas et all (|2011l) : 
and iDunklev et al.l (|2011f) . During roughly the same 
period, the South Pole Telescope also mapped the 
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microwave sky, and presented cosmolog i cal results in 
ILueker et all (l2^:IShlrokoff 'et all (l201lf):IKeisler et al.1 
(|201lD ; lReichardt et al.l (|2012f ): and lStorv et al.l (j2012ft . 

In this paper we describe a method to fit multi- 
frequency power spectra from the ACT data simulta- 
neously for CMB, foreground, and SZ p arameters, fol- 
lowing a si milar approach t o ana lyses in iDunklev et al.l 
(|2011l) and lReichardt et al.l (|2012D . We describe the like- 
lihood constructed for the 3-year ACT dataset, using 
data from the 2008-2010 observing seasons, and show 
how it can be used in combination with data from SPT 
in a self-consistent way. Using this likelihood from ACT, 
we then construct a simpler CMB-only likelihood, esti- 
mating CMB bandpowers marginalized over the SZ and 
foreground parameters. 

This is one of a set of papers on the 3-year ACT data; 
Das et ahl (2013) present the angular power spectra, and 
Sievers et al.1 ()2013h use the likelihoods presented here 
to estimate cosmological parameters. We begin in <J3] by 
describing the model for the mm-wave emission. In <|3] 
the full likelihood of the ACT data is described, including 
the combination with data from SPT and WMAP 7-year 
data. In Sj4] we show the small-scale model fit to the 
multi- frequency data. In <|5]we describe the compressed 
CMB-only likelihood, concluding in SjHl 

2. MODEL FOR THE MM-WAVE SKY 

Sky maps of mm-wave fluctuations at arcminute res- 
olution include components emitting at low redshift, in 
addition to the p rimordial CMB signal and secondary 
CMB effects (e.g..lSieyers et al.|[200a iFowler et al.| [201Ql: 
IDunklev et"atl 120111 : iKeisler et all l20l"Tl h" The power 
spectra from th e complete ACT dataset, reported in 
iDas et alJlpOll . are shown in Figure [IJ focusing on an- 
gular scales of interest for the primordial CMB signal. 
At scales smaller than a few arcminutes (£ ~ 1500) the 
secondary signal, which we define as the sum of fore- 
grounds and SZ effects, becomes significant compared to 
the CMB. We want to extract the primary CMB sig- 
nal, but, since there are more foreground components 
than frequency channels, information about both the fre- 
quency and scale dependence of the foregrounds is re- 
quired to separate the signals. In this section we de- 
scribe a model to fit the power spectrum of these fluctua- 
tions over the frequency range 90 v ^ 250 GHz probed 
by ACT, SPT, and other CMB experiments includin g 
the Planck satellite (jPlanck Collaboration I et al.| [2011). 
We follow a s i milar approach to iSievers et all (|2009f ); 
IDunklev et all pOlU ); iReichardt et all (|2012h . 

For frequency v and direction n we model the signal 
in the maps as 

ATO, n) = AT CMB (n) + AT SCC (V, n), (1) 

where AT CMB (n) are the lensed CMB fluctuations, 
which are independent of frequency in thermodynamic 
units. The secondary signal, AT sec (^, n), is dominated 
by the sum of tSZ and kSZ components, emission from 
dusty infrared galaxies and radio galaxies, and dust emis- 
sion from Galactic cirrus, all of which are functions of 
frequency. 

The cross-correlation power spectra between frequency 
Vi and Vj are calculated as 

cf = (f;( H )f e ( Vj )), (2) 



where Ti is the Fourier transform of T(n) in the flat-sky 
approximation. The theoretical cross-spectrum Bf 1 ' 11 
I : i + l)C* h ' ij /27r is modeled as 

where Bf MB is the lensed primary CMB power spectrum. 
In this analysis we model the secondary spectra as 

gsecij = £tSZ,ij + gkSZ.ij + gCIB-P.ij + gCIB-Cij 

+ £jtsz-ciB,ij + Sf d ' ij + S £ GaUj , (4) 

with contributions from the tSZ and kSZ effects; dusty 
galaxies that form part of the Cosmic Infrared Back- 
ground, both Poisson-like (CIB-P) and clustered (CIB- 
C); the cross-correlation between the tSZ effect and the 
CIB (tSZ-CIB); radio galaxies (rad); and dust emis- 
sion from Galactic cirrus (Gal). We assume that all 
other cross-spectra can be neglected. Measurements by 
WMAP and other CMB experiments show that other 
Galactic emission, including synchrotron and free-free 
emission, is negligible in the v £ 90 GHz freq uency range 
at th ese small scales and locations (e.g., iGold et all 
1 2 1 If ) . The cross-correlation of radio sources and both 
tSZ and CIB sour ces is also expected to be small (e.g., 
iSehgal et al.ll2010l who find a correlation of only a few 
per cent in simulations). Since the kSZ signal consists 
of positive and negative fluctuations, depending on the 
line-of-sight motion of the electrons that source the sig- 
nal, the two-point correlation function with other signals 
should average to zero. 

The majority of these secondary spectra will be com- 
mon to all regions of the sky and to different experi- 
ments. The power in the residual radio point sources is 
expected, however, to vary among data sets due to the 
removal of bright sources. For example, for a Poisson 
distribution of sources with differential number counts 
scaling as dN/dS oc S~ 2 , the Poisson power will be 

C e = J Q S ^S^dS, (5) 

i.e., Ci oc iSmax, where S max is the flux of the brightest 
sources in the map at a given frequency. Deeper surveys, 
with lower noise per pixel, are able to detect and mask 
out dimmer sources, so 5 max will be lower, leading to 
a lower residual power. Radio sources have a shallower 
dN/dS slope than CIB galaxies, so imposing a flux cut 
of, e.g., S'max = 15 mjy removes a significant amount of 
radio power, but little CIB power. 

The Galactic foreground power is also expected to vary 
between regions on the sky. The two regions mapped 
by ACT (Equatorial, ACT-E, and South, ACT-S) are 
shown in Figure [2] and summarized in Table [l] together 
with data from SP T. The temperatur e scale of the Galac- 
tic cirrus map from IFinkbeiner et al.l (|1999l ) is shown for 
comparison. A higher level of emission is expected in 
some regions of the ACT-E region. 

In the rest of this section we describe how each of the 
components in Eq. |4] are modeled. To allow for com- 
parisons between experiments, we normalize the power 
spectra at a pivot frequency of u Q = 150 GHz and scale 
£o = 3000. This differs slightly from t he convention used 
in previous analyses of the ACT data (jFowler et al.ll2010l : 
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Fig. 1. — Summary of small-scale mm-wave data measured by the Atacama Cosmology Telescope (Das ct al. 2013) and the South Pole 
Telescope (Kcisler ct al. 2011; Rcichardt ct al. 2012), in the angular range used for measuring the damping tail of the CMB. The ACT and 
SPT data are independently calibrated to WMAP. The vertical axis is i 4 C( instead of the con ventional £ 2 Ci to hi ghlight the features at 
these angular scales. The primary CMB signal corresponding to the best-fitting ACDM model (Sicvcrs ct al. 2013) is indicated (dashed), 
together with the total signal at 148 GHz (red, lower solid curve) and 217 GHz (black, upper solid curve), including secondary effects from 
SZ and foregrounds. Modeling the secondary contributions from SZ and foregrounds is vital to allow extraction of the primordial signal at 
small scales. 



TABLE 1 
Small-scale CMB datasets 



Dataset 


Frequency 3, 


Reference 


Area 


f ■ 


■^max 


& b 


ftsz c 


Z'Rad 


2^CIB 




GHz 




sq degrees 






mjy 


GHz 


GHz 


GHz 


ACT 


148 


Das et al. (2013) 


590 d 


500 


10000 


15 


146.9 


147.6 


149.7 




218 






1500 


10000 




220.2 


217.6 


219.6 


SPT-low 


150 


Kcisler ct al. (20111 


790 


650 


2000 


50 


152.9 


150.2 


153.8 


SPT-high 


95 


Reichardt et al. (2012) 


800 


2000 


9400 




97.6 


95.3 


97.9 




150 






2000 


9400 


6.4 


152.9 


150.2 


153.8 




220 






2000 


9400 




218.1 


214.1 


219.6 



a All cross-spectra between channels are used in the likelihood. 
b Flux cut imposed on map by its point sou rce mask. 

c Effective band-centers from ACT are from lSwetz et al.l <20TlD . given for tSZ, radio sources, and C1B sources. 
d This area includes the ACT-E region at dec = 0° (300 deg 2 ), and the ACT-S region at dec = -55° (290 deg 2 ). 



iDunklev et al.1 l20l"lT ) . In each case we describe the pa- 
rameterization used in the fiducial model; in a later sec- 
tion we consider possible extensions or modifications. 

2.1. Thermal Sunyaev-Zel'dovich 

Our model for the power from thermal SZ fluctuations 
is given by 

R tSZ,ij f{ v i)f{ v j) R tsz / fi \ 

B t - q tSZ J2(^ ) "0,1 i \ b ) 

where B^f is a template power spectrum correspond- 
ing to the predicted tSZ emission at vq for a model 
with amplitude of matter fluctuations erg = 0.8, nor- 
malized to 1 fiK 2 at Iq = 3000, and a t sz is a free pa- 
rameter describing its amplitude. An example is shown 



in Figure [H The factor f{v) = a;coth(x/2) — 4, for 
x = hv /ksTcMBi scales the expected tSZ emission to 
thermodynamic units at v, the effective band-center for 
the tSZ, gi ven in Table [H We ignore relativistic correc- 
tions fe.g-. lltoh et al.lll.998t) . since the low-mass clusters 
that dominate the spectra are well approximated by the 
non-relativistic for mula. This a t ,sz norm alization differs 
from that used in IDunklev et al.l (|2011[ ) . The present 
choice has the advantage of reducing the dependence on 
the choice of template, since the main difference between 
various templates is their amplitude. This means one 
expects to find the same constraint on atsz regardless of 
template, and a constraint on the SZ power can be con- 
verted back into a mo del- dependent constraint on erg. 
The template we adopt is derived from recent hydrody- 
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Fig. 2. — Regions of the sky used for ACT power spectra (dashed, 
IDas et al.|[20ll) in the Equatorial plane (ACT-E, 300 deg 2 ), and 
at -55° declination (A CT-S, 292 deg 2 ). The 800 deg 2 used for SP T 
power spectra (dotted, Keislcr et al. 2011; Rcichardt et al. 2012) is 
indicated, with 54 deg 2 o verlap with ACT-S. Th e color scales with 
Galactic cirrus intensity (Finkbciner et al. 1999). 

namic simulations described in Battaglia et al.l (|2012aj ). 
The simulations include the effects of radiative cool- 
ing, star formation, and feedback from AGN and su- 
pernovae. The predictions are consistent with SZ mea- 
surem ents from both SPT and ACT (e.g., iLueker et al.l 
2010), and the shape is shown in Figure [3] For the 
model with as = 0-8, th e predicted spectrum reported in 
IBattaglia et al.l (j2012al) has amplitude a t sz = 5.6 ± 0.9, 
with standard deviation estimated from ten simulations. 

Numerous other authors have also predicted the tSZ 
spectrum fro m independent simulations and analytical 
models (e.g. 
Seheal et al 



Komatsu fc Seliakl [20021: IShaw et al l 



20101: iTrac et al l 1201 It IShaw et al 



2009; 



2010; 



Battaglia et al] I2012ah lEfstathiou fc Migliacciol 120121) 
and the expected amplitude for fixed cosmological model 
varies depending on the astrophysical modeling of the 
clusters. However, the template shape is broadly consis- 
tent among models, and the data are not yet sensitive 
to shape difference, so we do not include a shape uncer- 
tainty. We do not mask clusters, and expect the total SZ 
power to be the same for ACT and SPT. 



Kinematic Sunyaev-Zel'dovich 
power is expected to have contribu- 



2.2. 
The kSZ 

tions arising from fluctuat ions in the electron density 
(jOstriker fc Vishniad 119861). and in the ionization frac- 
tion (e.g., iGruzinov fc Hul fl99l iMcQuinn etHI [20051 : 
llliev et al.ll2007f) . as well as from the motion of galaxy 
clusters at later times. We model the power as 



B 



kSZ,ij 



RkSZ 
flkSZ» £ i 



(7) 



where 23q^ z is a template spectrum for the predicted 
blackbody kSZ emission for a model with a a = 0.8, nor- 
malized to 1 fiK 2 at £o = 3000. The parameter aksz 
describes its normalization. We use a template that 
assumes a model with in s tanta neous reionization, de- 
scribed in Battagli a et al] (j2010) . This is derived from 
the same hydrodynamic simulations as the tSZ spectra 
in Sec 12. 1| and is shown in Figure [H The predicted am- 
plitude from the simulations is etksz = 1.5 for homoge- 
neous reionization at z = 10 in a a$ = 0.8 cosmology. 
This is a quarter of the expected tSZ power. The cor- 
re sponding kSZ template f or the 'nonthermal20' model 
in ITrac. Bode, fc Ostrikerl ((201 If ) has a similar ampli- 
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Fig. 3. — Template power spectra for t he thermal and kinetic 
Sunyaev-Zel'dovich effects (tSZ and kSZ, Battaglia et al. I2012al . 
20TcJl clustered CIB sources scaling as £ o s (CIB, lAddison et""aLl 
2012ri ). the cross-corr elation between tSZ and CIB (tSZ-CIB, 
neg ative at 150 GHz, lAdd ison et al. 2012c), and Galactic cir- 
rus ( Mivillc-Dcschencs & Lagachc 2005). They are normalized at 
I = 3000 and 150 GHz, and the tSZ-CIB is shown for a perfectly 
correlated signal. Poisson CIB and radio source power (not shown) 
scale as I 2 . 



tude and shape, as does the Sha w, Rudd. fc Na gai (2012) 
'CSF' model, and the IBode et al] ([20121 ) m odel. The 
power is expected to sc ale as roug hly of 5 " 5 (ITrac et all 
201l| [Shawet^[2012T) ). 



Reionization of the universe is not expected to be in - 
stantaneous, as was assumed by IBattaglia et all (2010). 
The shape and amplitude of the kSZ power from 
patchy reionization is far less certain, with simulations 
predicting a sign al at least as large as the homoge- 
neous signal (e.g.JZahn et aL1l2012l ; IMesinger et alj|2"012T : 



IBattaglia et alJl2012bD . This gives a total expected sig- 
nal of ~ 3 to 5 /iK 2 for simple reionization models at 
t = 3000, comparable to the tSZ at 150 GHz. The domi- 
nant effect of patchy reionization on the power spectrum 
at scales probed by ACT is to alter the amplitude, de- 
pending o n both the midpoint an d duration of reioniza- 
tion (e.g., IBattaglia et al.|[2012bl) . We test a modified 
shape in §4.3| but do not include additional shape uncer- 
tainty in the template in the basic model. 

2.3. Cosmic infrared background 

Thermal dust emission from high redshift star-forming 
galaxies, part of the Cosmic Infrared Background 
(CIB), is emitted in the rest-frame far i nfrared and 
redsh i fted into the mm- wave range (e.g., iPuget et al.l 
119961 : lHauser et al.l I1998D . Clustering of these galax- 
ies has been dete cted statistically in mm-wave maps at 
CMB frequency (l5airet~aTl [Mot iDunklev et all UOlTt 
Shirokoff et~aT1 120111 iPlanck Collaboration XVIII et al I 
20111 : lHajian et al.ll2012b iReichardt et aljl2012D. as well 
as in the sub- mm (e.g.. lLagache et al.l 120071: iViero et aTl 
2009D. Followin g the ana lyses in IDunklev et al.l ( l201ll)T 
Reichardt etall ([2012D ; lAddison et al] (|2012bD . the 



power from these galaxies is modeled as the sum of a 
Poisson and clustered component, given by 



B. 



CIB-P,ij 



i 



^(i/i,^ p )^(i/j-,/3p) 



/xK 2 (8) 



ACT Likelihood 
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for the Poisson part, and 



K,CIB-C.ij 

Be = a c 



to 



H{vi,f3 c )n(v 3l l3 c ) 
M 2 (^o,/3 c ) 



^K 2 (9) 



for the clustered part. Here, n is a power law index, and 
the frequency scaling of each component is given by a 
modified blackbody, 



f i(v,0) = t/ l B v (T d )g(v), 



(10) 



with emissivity indices /3 P and /3 C for the Poisson and clus- 
tered dust terms respectively. The function B u (Td) is the 
Planck function at frequency v for effective dust temper- 
ature Td, and the function g(y) = (dB u {T) / dT)^ 1 |t cmb 
converts from flux to thermodynamic units. The param- 
eters a p and a c normalize the two components at £q an d 
i/o, and different frequencies channels are assumed to be 
perfectly correlated. 

The frequency dependence we adopt is an approxima- 
tion to a sum of modified blackbodies at different red- 
shifts, so this emissivity and temperatu re are only ef- 
fective properties of the dust. Following lAddison et al.l 
(2012b]) we fix the dust effective temperature to Td = 
9.7 K. We also assume (3 P — (3 C in the basic model. 

The power-law angular scaling of the clustered term, 
with increasing £ 2 Ci power at small scales, is shown in 
Figure [31 and approximates the shape of the non-linear 
power spectrum, which includes contributions from pairs 
of galaxies in the same dar k matter halo, and bet ween 
galaxies in different halos. Ad dison et al.l ([2012b) find 
that a power law in £ provides a good fit to small- 
scale power spectra from Planck and the Balloon-borne 
Large- Aperture Submillimeter Telescope (BLAST), and 
from cross-correlating ACT and BLAST maps. This is 
consistent with observations of the correlation function 
from high-redshif t Lyman break galaxies, as well as lo- 
cal g alaxies (e.g., IGiavalisco et al.lll.998t IConnollv et al.l 
2002). We fix the power-law index to n — 1.2 in the 
fiducial case (^°' 8 ), i n close agreement with the estimate 
of n = 1.25 ±0.06 in lAddison et al1(|2012bfl . Both are in 
agreement with galaxy correlation functions. 

We do not expect the CIB power to vary significantly 
between the ACT and SPT maps, despite the different 
flu x cuts applied to rem ove sources. Using the model 
in lAddison et al.l ([2012aj ). the predicted effect of source 
masking on the CIB power is only at the per cent level. 

2.4. tSZ-CIB cross-correlation 

Some spatial correlation is expected between clusters 
that contribute to the tSZ, and CIB galaxies, since both 
trace the matter density field. The higher redshift and 
lower mass groups that make an important contribution 
to th e tSZ signal (lKo"matsu fc SeliaM I2002t iTrac et al.l 
1201 It [Battagli a et al. 2012bD are also likely to host dusty 
galaxies. lAddison et al.l ( 201 2d ) model this correlation, 
and predict the scale and frequency dependence of its 
angular power spectrum. For mm-wave spectra at £ > 
2000, a correlation of ~ 10 to 30% in power is predicted, 
with uncertainty dominated by uncertainties in the halo 
mass and redshift distribution of the CIB. A significant 
fraction of the CIB power on small scales is due to pairs 
of galaxies occupying group and cluster-mass halos (the 
'one- halo' term of the halo model). These same halos are 



responsible for the tSZ power. A tSZ-CIB correlation in 
units of power of tens of per cent is therefore possible on 
small angular scales even if the overall fraction of CIB 
emission associated with massive halos is small. 

The tSZ-CIB power is negative at 150 GHz, and can 
partially cancel power from the kSZ effect as it does not 
vary significantly with frequency over the range probed 
by ACT and SPT. As a result, neglecting this compo- 
nent ca n lead to artificia l ly tight constraints o n the kSZ 
po wer (IZahn et alJl2012t [Mesinger et al.M2012h . Follow- 
ing lAddison et al.l (|2012bf) we model the spectrum as 
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- 2/'K 



-B 



tSZ-CIB 
0.1 ' 



(11) 



where B e is the predicted correlation spectrum 

shape, normalized to 1 /iK 2 at £q and shown in Figure 
[3j The free parameter £ is the correlation coefficient. 
The Poiss on CIB parame t er a„ i s not included in E q. [Til 
unlike in iShirokoff etaD ([20111 ): IZahn et al.l (|2012f l. as 
the sources that dominate the CIB Poisson power in the 
mm-wave bands are unlikely to have significant red shift- 
overlap with the tSZ clusters ([Addison et al.ll2012d) . As- 
suming the same modified blackbody scaling for the CIB 
as in £12. 3i and tSZ frequency scaling f(v), the frequency 
scaling of the cross-spectra in thermodynamic units is 
then 



(12) 



for pivot scale v$. 

Since the correlation coefficient is poorly constrained 
by ACT, we impose a uniform prior of < £ < 0.2 in 
the basic model; the effect of widening the range to e.g., 
£ < 0.5, corresponding to the maximum allowed corr ela- 
tion in the models explore d by Addison et al. (20123), is 
discussed in iSievers et al.l ([20131 ). Due to the correlation 
between £ and the kSZ power, broadening the limit on 
£ increases the upper limit on the kSZ power, but does 
not affect cosmological results. 

2.5. Radio point sources 

The radio sources at ACT frequenci es are not expected 
to be significan tly clustered (see e.g., ISharp et alj|2010t 
lHall et aLll2010( ) , and to good approximation their power 
should be perfectly correlated between neig hbouring fre- 
quenc ies, consistent with simulations in Sehg aTet al.l 
(2010) and v alid for sources with the s ame spectral in- 
dices. As in iDunklev et al.l ([201 1[ ) and iReichardt et al.l 
(|20H . we model the residual power after masking bright 
sources as Poisson scale-free power, with 



B. 



rad,ij 



s 2 M 



/iK 2 (13) 



in thermodyamic units, where g{y) converts from flux 
units as for the CIB sources. The amplitude a s is nor- 
malized at vq and £q. Measurements of bright sources 
from ACT and SPT give an estimate for the spectral 
index in flux units of typi cally a s — —0.5 ([Vieira et al.l 
l2010t [Marriage et al. 20TTJ) . Assuming it holds at fainter 
fluxes, we fix a s = —0.5 in the fiducial model. 

Bright source counts can also be used to predict a s 
by extrapolating to fainter fluxes using a model for 
the number of sources as a function of flux. This was 
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do ne in iMarriage et al l (120111) for ACT and fo r SPT 
in iKeisler et all (|2Qlll ): iReichardt et ali (j2012D . Us- 



ing point so u rces m easured from the full ACT dataset, 
IGralla et all (pOll now predict a residual power a s — 
2.9 ± 0.4 after masking sources brighter than 15 mJy in 
both ACT regions (the level used to construct the mask 
for our maps in this analysis), where the catalog is es- 
timated to be complete. We impose this as a Gaussian 
prior on the power at 150 GHz. For comparison, the es- 
timated power in the SPT power spectra after masking 
to a flux level of 6 mJ y for SPT-high is a s = 1.3 ± 0.2 
(jReichardt et al.ll2012l ) , which we also impose as a Gaus- 
sian prior. 

2.6. Residual Galactic cirrus 

The G a lactic emission is spatially varying, and 
iDas et all (|2013l) show that dust emission can contribute 
significantly to the power spectr a, particularly in our 
Equatorial region. As reported in IDas et all (|2013l ). we 
apply a mask to regions of high dust emission before 
computing the pow er spectrum, using measurements a t 
100 urn from IRIS (jMiville-Deschenes fc Lagachell2005D . 

We then marginalize over a residual Galactic cirrus 
component using a power-law template 



B 



Gal.ij 



tig ■ 



frequency index j3 g 



yK\ (14) 



and 



-0.7. This angular scal- 



with amplitude 
angular scaling 
ing is estimated from the 100 nm IRIS dust maps 
(jMiville-Deschenes" fc Lagachc 2005). The frequency 
scaling is e stimated by cor relating the IRIS dust maps 
with ACT (|Das et aX1l2013D, and is consistent with early 
resul ts from Planck ( Planck C ollaboration XIX et all 
20111). U s ing th e correlation coefficients estimated in 
Das et all (|2013l ). we impose priors of a ge = 0.8 ± 0.2, 
and a gB = 0.4 ± 0.2 in the ACT-E and ACT-S spectra 
respectively. 

For the SPT data, a small Galactic cirrus residual is 
also exp ected. In our b a sic m odel we follow the treat- 
ment in IReichardt et all ([2012D . fixing the Galactic cir- 
rus power to S3000 = 0.16, 0.21, and 2.19 yK 2 in ther- 
modynamic units at 95, 150, and 220 GHz, with scale 
dependence cx I -1 - 2 . However, this model has a steeper 
angular power law than in our ACT model, and a shal- 
lower frequency scaling ((3 = 3.6 between 150-220 GHz). 
For consistency we therefore test the effect of adopting 
the ACT model instead, using £~ - 7 and /3 g = 3.8, and a 
prior of a g = 0.4 ± 0.2. We find no effect on parameters 
and no change in the goodness of fit. 

3. FULL LIKELIHOOD FROM SMALL-SCALE DATA 

In this section we describe the multi-frequency likeli- 
hood used to model the ACT data, and show how we 
extend it to include other small-scale datasets, in partic- 
ular data from SPT. 



fro m ACT maps obtain ed using the method described 
in iDunner et all ([2012). The likelihood is Gaussian- 
distributed to good approximation. To construct the 
likelihood for each region, given some model spectra 



, we compute bandpower theoretical spectra using 
= wliCf 1 ' 13 , where w l A is the bandpower window 



function in ban d b for cross-spectrum ij, described in 
IDas et all (poll . 

The likelihood, Jz? , of the data for each ACT region 
separately is given by 



21nJ^= (C, 



-C b )+ In detS, (15) 



where E is the bandpower covariance matrix. Each of the 
model and data vectors Cj h and C b contain three sets of 
spectra, 

n r^l48,148 z-,148,218 ,-,218,2181 

= [C 6 ,C b ,C b j, (16) 

for ACT-E and ACT-S separately, and each spectra set 
C % b 3 itself contains spectra for each cross-season. There 
are two seasons used for ACT-E (3 cross-season spectra), 
and three for ACT-S (6 cross-season spectra) . The total 
likelihood is given by 

- 21ni?ACT = -21n^ AC T-E - 21nJSf A CT-S- (17) 

3.1.1. Calibration and beam uncertainty 

The data power spectra are calibrated, but have un- 
certainties. We therefore include a calibration parameter 
Ui, for each map i, that scales the estimated data power 
spectra as 

-> ViVjC'l (18) 
and the elements of the bandpower covariance matrix as 



^bb' 



(19) 



To account for both ACT regions, we include four cali- 
bration parameters: t/i e , y2e for ACT-E at 148 GHz and 
2 18 GHz, and y^, y 7fl for ACT-S. 

IDas et all (|2013f ) calibrate the 148 GHz maps using 
WMAP, following the method in lHaiian et all (|201lD , 
at an effective £ = 700, resulting in a 2% map calibration 
error in CMB temperature units. We impose this as a 
Gaussian prior, with yi e ,y u = 1.00 ±0.02. The 218 GHz 
maps are calibrated relative to 148 GHz, at an effective 
£ = 1500. The 218 GHz calibration is constrained by the 
cross-spectrum, so no prior is imposed on yi e and yis- 
Within each frequency, the individual seasons are cali- 
brated to each other; the inter-season calibration error is 
absorbed into the single overall calibration uncertainty. 

Uncertainties in the measured beam window functions 
for ACT at 148 GHz are between 0.7 and 0.4%, and at 
218GHz between 1.5 and 0.7%. We incorporate uncer- 
tainties in the measured beams by including them di- 
re ctly in the c ovaria nce matrix for the spectra, described 
in IDas et all ((2013). This technique assumes a fiducial 
model for the power spectra but is insensitive to its exact 
form. 



3.1. Likelihood from the ACT data 

The data from IDas et all (|2013l ) describe the two 
ACT regions separately (ACT-E and ACT-S); and con- 
sist of multi-season and multi-frequency spectra, with 
an associated covariance matrix. They are derived 



3.1.2. Secondary model parameters 

Our model described in Section [5] has nine free sec- 
ondary parameters for ACT in the basic case: a t sz and 
flksz describing the SZ emission, a p , a c and j3 c describ- 
ing the CIB power, a s describing the radio power, £ 
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nent; we use the band-c enters for SZ, radi o, and dusty 
sources given in Table 1 (ISwetz et al.H201lD . 



Fig. 4. — (Top) Power spectra measured by ACT l|pas ct al. 20TJ) at 148 and 218 GHz, and their cross-spectrum, coadded over ACT-E 
and ACT-S. We show the primary (lensed CMB in dotted black line) and secondary contributions (dotted lines) to the best-fitting model. 
(Bottom) Residual power in the ACT cross- frequency spectra, after subtracting the best-fitting model, at 148 (left), 148x218 (center), and 
218 GHz (right). The errors at small scales are correlated due to beam uncertainty. The model is a good fit simultaneously to ACT-E and 
ACT-S, with no sigificant residual features. 

describing the tSZ-CIB cross-correlation, and a ge and 
a gs describing the Galactic cirrus emission. The lat- 
ter four have strong priors imposed, as described in <j2j 
a s = 2.9 ± 0.4, < £ < 0.2, a ge = 0.8 ± 0.2, and 
a gs = 0.4 ± 0.2. In addition to the nine model param- 
eters, there are four calibration parameters for ACT. In 
£14.31 we investigate how additional, or fewer, parameters 
affect the fit of the model to the data. To compute the 
model requires an effective frequency for each compo- 



3.2. Combining with SPT data 

The South Pole Telescope ob served the sky f rom 2 007- 
10. Spectra are reported in iKeisler et al.l (|201 lh for 
angular scales 650 < I < 3000 at 150 GHz, and in 
iReichardt et al.l (|2012D for angular scales 2000 < I < 
9400 at 95, 150 and 220 GHz. These observations are 
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Fig. 5. — Distributions for secondary parameters from ACT and SPT, for best-fitting ACDM model. Parameters {atsz, <2kSZ> a p> a c, 
a gs , a ge , as} are the £(£ + 1)C(/2tt power in /^K 2 at t = 3000 and frequency 150 GHz. The tSZ-CIB correlation parameter £ is also defined 
at i = 3000. The dust emissivity index /3 C is in flux units, for a mod ified b lackbody with effective temperature 9.7 K. Conversions to power 
at each frequency are given in Table[2] Strong priors, described in i|3.1.2l are imposed on {£, a ge , a gs , a s }. 



summarized in Table [T] One of the goals of our work 
here is to test for consistency between the two exper- 
iments, by using a common framework to describe the 
SZ and foreground components. As this article was be- 
ing prepare d, refined sp e ctra f rom SPT at 150 GHz were 
reported in iStorv et al.l ([2012D ; we do not include these 
latest data in our comparison. 

Before fitting the SPT data with the ACT secondary 
model, we confirm that we reco ver the parameters and 
X 2 obtained for the model used in lReichardt et all (|2012f ) 
using the SPT data. To combine the d ata over the full 
angula r range, we f ollow the met h od in |Rcicliard t~et al.l 
(|2012ft , using the iKeisler et all (1201 ID data at i < 
2000 (SPT-low) and the IReichardt et all (|2012ft data at 
smaller scales (SPT-high). More radio source power has 
been removed from the SPT-high spectra due to mask- 
ing at a deeper flux level, so the expected residual radio 
power in SPT-low is 63000 = 10.5±2.4 /iK 2 , compared to 
1.3 ± 0.2 for SPT-high. We account for this by first sub- 
tracting a radio Poisson power of Bg = 9 .2 /iK 2 from the 



SPT-l ow data, following the approach in lReichardt et al.1 
(2012). A Gaussian prior is then imposed on the overall 
residual radio level in SPT of a s > — 1.3 ± 0.2. 

We then extend the ACT secondary model to fit the 
SPT power spectra. Six of the ACT model parameters 



are expected to be common for the SPT data (the SZ 
and CIB parameters: chsz, Oksz, £, o, p , a c and j5 c ). In 
addition, to fit the SPT data we require a separate radio 
source parameter, a s > , and three calibration parameters, 
2/2, 2/3, to calibrate the 95, 150, and 220 GHz maps 
respectively. We impose a uniform prior on these calibra- 
tion parameters, as the SPT covariance matrices include 
the calibration uncertainty. 
The likelihood for ACT and SPT together is given by 



2 In JSf = -21nif AC T - 21ni? S PT 



(20) 



The SPT likelihood is constructed as in Eq. [15J with 
model and data vectors 



C b 



r.^,95,95 ^95,150 ^95,220 .-,150,150 ^150,220 ^220.220] 
l U h >°b '°1 i°>. 



(21) 

for SPT-high at I > 2000, and C b = C* fc 150 ' 150 for SPT-low 
at t < 2000. To compute the model, we use the band- 
cent ers for SZ, radio and du sty sources given in Table 1, 
from IReichardt et al.l (|2012ft . 

There is some degree of covariance between the ACT 
and SPT spectra, due to the 54 deg 2 overlapping region 
of sky. Covariance between the two spectra due to cos- 
mic variance, scaling as l// s ky, is estimated at a level of 
8% (using 54/V590 x 800); the addition of noise lowers 
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this level, so we neglect the correlation in our combined 
analysis. 

3.3. Multi-frequency likelihood prescription 

To return the ACT (or ACT+SPT) multi-frequency 
likelihood for a given model we follow this approach: 

• Select primary cosmological parameters, and com- 
pute a theoretical lenscd CMB power spectrum 
gCMB ugmg the C AMB numerical Boltzmann code 
(|Lewis et al.|[2000l) . 

• Select values for common secondary parameters: 
= {atsz, OkSZ, a p , a c , Pc\. 

• Select values for ACT-specific secondary and cal- 
ibration parameters: 9 = {a s , a ge , a gs , yi e , ?/2e, 
Vis, J/2s}, and/or SPT-specific secondary and cali- 
bration parameters: 9 — {a s >, yx, y2, 2/3}. 

• Compute the total theoretical secondary power 
spectra £^ CC,IJ for all the required cross-spectra with 
Eq. |4j using the effective frequencies for each exper- 
iment. 

• Compute the total model power at each frequency, 

gth.ij = g CMB + B |ec,ij_ 

• Compute the bandpower theoretical power spectra 
for each dataset for both South and Equatorial re- 
gions for ACT (and for SPT), and compute the 
likelihood using Eq. [TTl 

3.3.1. Combining with WMAP 

In lSievers et~aT1 (|2013ft we use the ACT and SPT like- 
lihood in combination with data from WMAP seven-year 
data to estimate cosmological parameters. The WMAP 
data measure t < 1000 angular scales, and so have min- 
imal contamination from SZ and point sources. The 
public 7-year likelihood estimates the temp erature spec- 
trum from V and W bands (61 and 94 GHz. lLarson et all 
120111 ). At these frequencies and angular scales, the in- 
frared point source contribution is expected to be negli- 
gible, consistent with ACT and SPT measurements. The 
radio point source level is estimated and subtracted inter- 
nally to th e WMAP analysis using the multi-frequency 
data (e.g., iNolta et al.ll2009D . Finally, in light of ob- 
servations by both ACT and SPT, we also neglect the 
SZ power in the WMAP data, as it is expected to be 
small: 1(1 + DC, nnn/2? r ~ 12 \x¥? at 61 GHz from the 
iBattaglia et~atl (|2012aD model, assuming as — 0.8. 

4. TESTS OF THE MULTI-FREQUENCY LIKELIHOOD 

In this section we test the goodness of fit of the model 
to the ACT power spectra, assuming the ACDM cosmo- 
logical model. We estimate the probability distributions 
of the secon dary paramete r s usin g the MCMC method 
described in lDunklev et all (|2011[ ). fixing the ACDM pa- 
rameters at best-fitting valuefl We then investigate a 

1 Physical baryon density VL^h? = 0.02226, cold dark matter 
density Q c h 2 = 0.1122, ratio of the acoustic horizon to the angu- 
lar diameter distance at decoupling = 1.040, scalar amplitude 
ln[10 10 A s ] = 3.186, spectral index n a = 0.9707, and optical depth 
r = 0.898. 



set of possible extensions or modifications to the sec- 
ondary model. We include the SPT power spectra and 
examine the consistency of the foreground model between 
the two datasets. 

4.1. ACT data 

We find that the model provides a good fit to the ACT 
data over the full range of angular scales and frequen- 
cies. Figure 0] shows the total spectra (coadded over 
ACT-E and ACT-S, with the best-fitting Galactic cirrus 
component removed) decomposed into primary and sec- 
ondary contributions. The SZ and foregrounds dominate 
at I > 2400 at 218 GHz, and at I > 3200 at 148 GHz. 
The goodness of fit is \ 2 = 675 for 697 dof (reduced 
X 2 = 0.98, with PTE= 0.72, for 710 data points and 13 
parameters) . Figure [4] also shows the residual power af- 
ter subtracting the best-fitting model; we do not observe 
any significant features, indicating that the model fits 
both the angular and frequency dependence of the data 
in both regions. 

The marginalized distributions for the secondary pa- 
rameters fitting the data are shown in Figure [5] and sum- 
marized in Table [2] The Poisson-like and clustered CIB 
power, a p and a c , are detected at high significance, with 
index /? = 2.2±0.1 consistent with lAddison et all (|2012bl ) 
who find 2.20 ± 0.07. The tSZ and kSZ power are indi- 
vidually seen at low significance, with an anti-correlation 
between atsz and aksz- The kSZ power peaks at a non- 
zero value, but the distribution is broad and consistent 
with zero. The total SZ power is detected at high sig- 
nificance. The tSZ-CIB correlation coefficient is uncon- 
strained in the prior range < £ < 0.2, and is also uncon- 
strained by ACT if allowed to vary over a broader range 
(e.g., £ < 0.5). The parameters for the power from radio 
sources and fr om Galactic c i rrus a re driven by their prior 
distributions. iSievers et al.l (2013) present a physical in- 
terpretation of these parameters and their degeneracies; 
the constraints are c onsistent with those f ound in the 1- 
year ACT analysis in lDunklev et al.l ()201lD . with reduced 
errors. 

In Figure [6] we show the individual components that 
contribute to the 148 GHz and 218 GHz power spectra 
after removal of the best-fitting Galactic cirrus power. At 
148 GHz there are contributions from all the components. 
At 218 GHz the secondary spectrum is dominated by 
dusty point sources, both clustered and Poisson. This is 
illustrated further in Figure[7J which shows the frequency 
dependence of the dominant components in our model at 
£ = 3000. The derived constraints on the CIB and radio 
source components, and the Galactic cirrus emission, at 
the ACT effective frequencies are also given in Table [3] 
to allow comparison with other models. 

4.2. Combination with SPT 

The same model also provides a good fit to the SPT 
spectra. The SPT data extend the frequency range to 
95 GHz, adding three additional cross-spectra to the like- 
lihood. We show the parameters estimated from SPT 
alone in Figure [5] and in Table [2j they are consistent with 
those from ACT, with ~ la shifts in aksz and /3 C . The 
radio Poisson level is lower due to the greater number 
of radio sources masked in the SPT maps. The model 
is shown with the SPT spectra in Figure [8j the good- 
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Fig. 6. — Power spectra measured by ACT at 148-218 GHz, with the best-fitting individual SZ and foreground components from Table 
[2] The Galactic cirrus component has been subtracted. At 148 GHz (top) the secondary components are significant at scales smaller than 
I ~ 2000, with contributions from tSZ, kSZ, radio galaxies, the CIB, and the tSZ-CIB cross-correlation. The tSZ, kSZ and tSZ-CIB are 
non-zero in this model, but are not individually significantly detected from the ACT spectra. The radio power is constrained by bright 
source counts. At 218 GHz (bottom) the secondary signal is significant by i ~ 1000, and is dominated by the Poisson and clustered CIB. 



ness of fit is X 2 = 96 for 107 dof (reduced x 2 = 0.89, 
PTE=0.77). 

Given the consistency of the two datasets, we combine 
them to generate a joint likelihood; Figure [5] includes 
the secondary parameters derived from a joint fit. In 
this case there are ten foreground parameters, and seven 
calibration parameters. The tenth foreground parame- 
ter (not plotted) is a s i for the Poisson radio sources in 
SPT. The goodness of fit of the joint model is \ 2 = 773, 



which can be compared to \ 2 = 675 + 96 = 771 for the 
independent fit to each data set. This supports their 
consistency. 

We report the derived constraints on the CIB and ra- 
dio source components at the ACT and SPT effective 
frequencies for each band in Table [3] A difference of 
approximately 15% is expected between the CIB power 
at 148 GHz for ACT and 150 GHz for SPT, due to dif- 
ferent effective bandpass frequencies and the strong CIB 
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TABLE 2 

Likelihood parameters, assuming best-fit 6-parameter ACDM for the lensed CMB a 





Parameter 


Prior b 


ACT C 


SPT 


ACT+SPT 


sz 


atsz 


> 


3.3 ± 1.4 


4.1 ± 0.9 


4.0 ± 0.9 




a kSZ 


> 


< 8.6 


< 4.2 


< 5.0 d 


CIB 


l-lp 


> 


6.9 ± 0.4 


7.0 ± 0.4 


7.0 ± 0.3 




a c 




A Q 4- Q 




k 7 -t- n fi 




Pc 


> 


2.2 ±0.1 


2.0 ±0.1 


2.10 ±0.07 


tSZ-CIB 


e 


< i < 0.2 


< 0.2 


< 0.2 


< 0.2 


Radio 




2.9 ±0.4 


3.1 ±0.4 




3.2 ±0.3 




a s i 


1.3 ±0.2 




1.4 ±0.1 


1.4 ±0.1 


Galactic cirrus c 


Uge 


0.8 ±0.2 


0.9 ±0.2 




0.9 ±0.2 




dgs 


0.4 ±0.2 


< 0.73 




< 0.70 



Calibration y le 1.00 ±0.02 1.010 ±0.007 1.006 ±0.006 

y 2e 0.99 ±0.01 0.99 ±0.01 

y u 1.00 ±0.02 1.011 ±0.007 1.010 ±0.007 
y 2s 1.03 ±0.01 1.02 ±0.01 

yi 1.01 ±0.02 1.01 ±0.02 

y 2 1.007 ±0.008 1.008 ±0.008 

2/3 — — 1.02 ±0.02 1.03 ±0.02 

best fit x 2 /dof 675/697 96/107 773/810 

PTE 0.72 0.77 0.82 



a Secondary parameters marginalized over the 6 ACDM model parameters are reported in Table 1 offSievcrs et al. (2013), and are consistent 
with these results. The marginalization has little effect on these secondary parameters, increasing errors by at most 10%. 
b A flat prior is imposed, unless indicated as a Gaussian with x ± y for mean x and standard deviation y. 
c R,esults are reported as 68% confidence levels or 95% upper limits; £ is unconstrained so the prior upper limit is reported. 
d If the prior on § is broadened to < £ < 0.5, the upper limit increases to a^sz < 6.9 {Sicvcrs ct al. 20l3). 

c The SPT cirrus level we use is B s000 = 0.16, 0.21, and 2.19 /xK 2 at 95, 150, and 220 GHz, as measured in IReichardt et al.l (I20T2T) . 



TABLE 3 

Derived constraints on foreground power, H3000 (a*K 2 ) 



ACT 

148 GHz 218 GHz 


95 GHz 


SPT 
150 GHz 


220GHz 


CIB-P 
CIB-C 


6.8 ±0.4 
4.8 ±0.9 


78 ±12 
54 ±16 


0.90 ±0.02 
0.76 ±0.02 


8.0 ±0.5 
6.8 ±0.8 


69 ± 10 
59 ± 12 


Radio 


3.2 ±0.4 


1.4 ±0.2 


7.2 ±0.8 


1.4 ±0.2 


0.7 ± 0.1 


Gal-E a 
Gal-S 


0.9 ±0.2 
0.4 ±0.2 


11 ± 2.3 
5.0 ±2.3 









a Gal-E and Gal-S are the Galactic cirrus powers in the ACT-E 
and ACT-S spectra. The levels are close to t he priors impose d 
from the measured cross- correlations with IRIS HDas et al.ll2013ft 

frequency dependence across the mm- wave bands. 

4.3. Tests of the likelihood 

This model fits the ACT and SPT data, and includes 
our uncertainties about the physical components, with 
priors describing our knowledge from other observations. 
However, it is a simplified parameterization of the emis- 
sion. We therefore consider a set of extensions or modi- 
fications to the model, and test how the goodness of fit 
to the ACT data is affected by an increase or decrease 
in parameters, or a change in the prior assumptions. In 
these tests, summarized in Table 21 we hold the cosmo- 
logical model fixed at the best-fitting ACDM parameters. 
A subset of t hese extensions are considered further in 
ISievers et al.l (|2013[ ). testing their effect on the primary 
cosmological parameters. 




140 160 180 
v [GHz] 



Fig. 7. — Frequency dependence of the dominant components 
of the foreground power at I = 3000 measured by the combined 
ACT and SPT data sets. The bands show the 1<t uncertainties 
from Table [2] At 150 — 220 GHz the power from fluctuations in 
the CIB dominates; at lower frequencies the thermal SZ and radio 
source power is more significant. The SPT radio power is lower 
due to deeper integration. The kSZ and tSZ-CIB components are 
not shown. 

The CIB appears to be well-fit current ly by a power- 
law in angular scale, with Bi oc £ 8 . lAddison et al.l 
(20121)) find an uncertainty of 0.0G in this scaling. If we 
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Fig. 8. — Power spectra at 95 , 150 and 220 GHz, and their 
cross-spectra, measured by SPT (Rcichardt ct al. 2012), fit with 
the same model as the ACT data in Figure [4] At 150 GHz the 
SPT-low spectrum from Kcislcr ct al. (2011) is included, with ex- 
cess radio power subtracted for comparison. Accounting for the 
different flux cuts applied to the ACT and SPT maps, and the 
different bandpass effective frequencies, the spectra are consistent. 



TABLE 4 
Modifications to secondary model 



Model 


Number of 


ACT b 




parameters 3. 


A X 2 


Fiducial 


9 





CIB index n free 


10 


-4 




10 





CIB Poisson corr = 0.8 


9 


5 


CIB Clustered corr = 0.8 


9 


1 


CIB T d = 13.6 K 


9 


1 


Fixed kSZ, a kSZ = 1.5 


8 


3 


Altered kSZ shape 


9 


1 


No tSZ-CIB corr, £ = 


8 


2 


No SZ 


6 


21 


Radio index a a = 


9 


1 


No Galactic residual 


7 


6 



a Does not include calibration parameters. 
b We report the A% 2 to the nearest integer. 



allow the index n to vary, we find no improvement in the 
fit, but parameter distributions for the CIB parameters 
dp and a c are broadened, as they are correlated with the 
power-law scaling. We also assume that the CIB emission 
is perfectly correlated among frequencies, in the range 
95-220 G Hz. Evidence for imperfect correlation w as re- 
ported in IHIncjn^oiia^ The 
effect of this assumption is tested by setting the corre- 
lation coefficient to < 1 in the model, choosing 0.8 for 
either the Poisson or clustered components, roughly cor- 
respondi ng to the degree of correlation between ma ps re- 
ported in lPlanck Collaboration XVIII et all (|2011[ ). For 
the Poisson component we find that this degrades the 
goodness of fit by A\ 2 = 5 compared to the perfectly 
correlated case. We also assume a common frequency 
scaling of the clustered and Poisson terms. A different 
scaling in frequency may be expected if, for instance, the 
redshift dependence of the clustered and Poisson power is 
different, so a common index is not necessarily expected. 
Allowing it to vary independently does not significantly 
improve the goodness of fit, but does lead to a poorly 
constrained distribution for /3 C , and removes the detec- 
tion of the CIB at 148 GHz. The index for the Poisson 
sources in this case, (3 p = 2.14 ± 0.15, is consistent with 
the joint index. Changing the effective dust temperature 
fro m 9.7 K to 14 K, co nsistent with the value obtained 
bv lGispert et~aLl (|2000l ) from a fit to the FIRAS CIB fre- 
quency spectrum, has no effect on the model, apart from 
a corresponding change in j3 c . 

Exploring the SZ assumptions, we consider fixing the 
kSZ contribution to the linear theory estimate for a uni- 
verse with as — 0.8, assuming homogeneous reioniza- 
tion. This degrades the goodness of fit of the model 
by only A% 2 = 3, indicating that the data cannot yet 
distinguish between homogenous and patchy reioniza- 
tion. Limiting the kSZ in this way also leads to tighter 
constraints on the tSZ power. Modifying instead the 
shap e, we find that adding a p atchy reionization template 
from lBattagli a~et al.l (|2012bl ) to the kSZ, which changes 
its shape, does not affect the other secondary parame- 
ters. Neglecting the tSZ-CIB cross-correlation also fits 
the data equally well for one fewer parameter. The de- 
pendence of the kSZ constraints on the model for the 
tSZ-C IB correlation is explored further in iSievers et al.l 
(|2013|) . If we neglect the SZ components altogether, set- 
ting atsz = »ksz = and keeping only the CIB and 
Galactic components, the goodness of fit significantly 
worsens, with an increase of A\ 2 — 21. 

Our model imposes an a priori assumption on the fre- 
quency scaling of the extragalactic radio sources. We test 
the effect of changing the radio spectral index to a s =0, 
finding negligible effect on parameters and goodness of 
fit. There is al so little effect from changing the prior on 
the power from lGralla et al.l (|2013t ) by la, corresponding 
to a different model for the bright radio sources that lie 
above the detection threshold. Removing the prior alto- 
gether opens up degeneracies with other parameters, but 
does not significantly improve the goodness of fit. 

Finally, we test the effect of removing the Galactic cir- 
rus components; the goodness of fit worsens by A% 2 = 6 
and the clustered CIB level increases. This indicates a 
preference for Galactic cirrus at the 95% confidence level. 



5. CMB-ONLY LIKELIHOOD 
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Fig. 9. — Estimated CMB bandpowers from ACT, marginalized over extragalactic source and SZ components. Bandpowers are estimated 
for ACT-E and ACT-S separately; here we show the inverse-variance weighted combination. The bandpowers are correlated at the ~ 20% 
level at scales I ^1 2000 due to covariance with the secondary parameters. The total multi-frequency spectra for ACT-E (dashed, at 148 GHz, 
148x218 GHz, and 218 GHz) are also shown to indicate the significant level of SZ and foreground power at small scales. 

5.1. Method: bandpowers via Gibbs sampling 

To implement this method in practice, we estimate rib 
CMB bandpowers, marginalizing over the secondary pa- 
rameters. We use the full multi-frequency likelihood from 
Sj3]to do this, but estimate CMB bandpowers instead of 
cosmological parameters. 

We recall that the model for the theoretical power for 
a single cross- frequency, cross-season spectrum, C^ h:I \ is 
written as 



Understanding the contribution of the secondary com- 
ponents to the ACT power spectra is vital for extracting 
the cosmological information, due to possible degenera- 
cies between primary and secondary parameters. Values 
of the secondary parameters are also astrophysically in- 
teresting. However, if we are only interested in the cos- 
mological parameters, a simplified likelihood is desirable. 

We construct a CMB-only likelihood from ACT data 
as follows. Instead of using the ACT likelihood to es- 
timate cosmological parameters, we take the interme- 
diate step of estimating the CMB power spectrum in 
bandpowers, marginalizing over the possible contamina- 
tion. This is a natural extension to forms of CMB data 
compress ion that have been adopte d in earlier analy- 
ses (e.g.. iBond. Jaffe. fc KnoxH ^OOO). Such 'grand uni- 
fied spectra' were used in a number of subsequent pa- 
pers to combine the results from various CMB experi- 
ments, margin alizing over a variet y of nuisance param - 
eters, e.g. in iSievers et al.1 (|2003h : IBond et al.1 (|2003| ): 
IBond et all (120041) . At large scales, where contamination 
from SZ and point sources is negligible, the estimated 
CMB is simply an optimally co mbined average of th e 
multi- frequency spectra as in e.g., Iffinshaw etafl (|2003l) . 
At smaller scales the CMB spectrum has additional un- 
certainty due to secondary contamination. 

By marginalizing over nuisance parameters in the 
spectrum-estimation step, we can effectively decouple the 
primary CMB from non-CMB information. No addi- 
tional nuisance parameters are then needed when esti- 
mating cosmological parameters. 



C 



th.ij _ pCMB 



CSCC.1J 
£ 



(0), 



(22) 



where C| ec ' 1J (#) is the secondary signal as in Eq. 21 and 
is a function of secondary parameters 9. Writing the 



spectrum in bandpowers, C t 



th.ij _ 



ffeCl ' l K where 



w are 



the bandpower window functions, we write the model for 
the bandpowers in vector form as 



Cf = ACF MB 



+ cr(0), 



(23) 



where C^ h and C| cc are multi-frequency, multi-season 
vectors of length rib x n spe c, where rib is the number of 
bandpowers, and n spoc is the number of cross-season and 
cross-frequency spectra (n spcc = 9 for ACT-E, and 18 
for ACT-S). The secondary spectra differ between fre- 
quencies but not between seasons. The mapping matrix 
A, with elements that are either 1 or 0, maps the CMB 
bandpower vector (of length n&), which is the same at 
all frequencies and in all seasons, onto the (rib x ^spec) - 
length data vector. 

We want to estimate C^ MB , marginalized over the sec- 
ondary parameters, 9. The posterior distribution for 
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MB j given the observed multi- frequency, multi-season 
spectra C b , can be written as 

p{C™ B \C b ) = Jp(C™ B ,9\C b )p(9)d9. (24) 

We find that Gibbs sampling provides an efficient way to 
map out the joint distribution p(C b MB , 0\Cb), and to ex- 
tract the desired marginalized distribution p(C b MB |Cf,). 

Gibbs sampling can be used in the special case that at 
least one conditional slice through a multi-dimensional 
distribution has a known form, and has been used, for 
example, to estimate the large-scale CMB power spec- 
trum, and to margina l ize over Gala c tic fo r egrounds (e.g., 
Wandelt et all 12004 Uewell et all 120041: Eriksen et al.1 
2004HDunklev et al.ll2009HLarson et alJl2011h . Here, we 
split the joint distribution into two conditional distribu- 
tions: p{C^ MB \9,C b \ and p{9 | C fc CMB ,C b ). We write the 
multi-frequency likelihood for a single ACT region, from 
Eq. [H as 

- 2 InJz? = (AC b CMB + C b scc - C 6 ) T S- 1 (AC b CMB + C s b cc - 

+ lndet£, 

which is a multivariate Gaussian. 

If C| cc is held fixed, the conditional distribution for the 
CMB bandpowers, p(C b MB \9,C b ), assuming a uniform 
prior for p(C b CMB ), is then also a Gaussian. It has a 
distribution given by 

- 2\np(C? MB \e,C b ) = (C™ B - C b fQ-\C™ B ~ C b ) 

+ lndetQ, (26) 

The mean, C bl and covariance, Q, of this conditional 
distribution are obtained by taking the derivatives of the 
likelihood in Eq. [55] with respect to C b CMB . This gives 
mean 

C b = [A T S- 1 A]- 1 [A T S- 1 (C b - CH], (27) 
and covariance 

Q = A T £ *A. (28) 

We can draw a random sample from this Gaussian dis- 
tribution by taking the Cholesky decomposition of the 
covariance matrix, Q = LL T , and drawing a vector of 
Gaussian random variates G. The sample is then given 
by C b CMB = C b +L- 1 G. 

If instead C b MB is held fixed, the conditional distri- 
bution for the secondary parameters, p(9\C b MB , C b ), is 
not a Gaussian, but can be sampled with the Metropo- 
lis algorithm that is used in the MCMC sampling in <21 
To map out the full joint distribution for 9 and C b MB 
we alternate a Gibbs sampling step, drawing a new vec- 
tor of CMB bandpowers, C b , with a Metropolis step, 
drawing a trial vector of the secondary parameters 9. 

We choose a uniform positive prior distribution for 
p(C b CMB ), and restrict the CMB bandpowers to be zero 
at £ > 4500, where the CMB power is expected to 
be less than 1 /xK 2 . About 100,000 steps are required 
for convergence of the jo int distribution, assessed with 
the IDunklev etlH (<2005D spectral test. The mean and 
covariance of the resulting marginalized bandpowers, 
C b CMB , are then esti mated following the st a ndard MCMC 
presc ription (e.g., iLewis fc Bridle) [2001 ISpereel et all 
12001 . 



5.1.1. Combining spectra from different regions 

There is only one underlying CMB power spectrum, so 
this method could be used to estimate a single spectrum, 
or set of bandpowers, from the two ACT regions. How- 
ever, the bandpower window functions are different for 
each region due to their distinct geometries. To easily 
conserve this information, we estimate the CMB band- 
powers for ACT-E and ACT-S separately. Since the sec- 
ondary parameters are common to both, the estimated 
CMB bandpowers will be correlated between the regions 
at small scales. 

To estimate the joint distribution for the 
ACT-E and ACT-S bandpowers, we map out 
p{Cf MB - E ,Cf MB - s ,9\C b ) by taking sequential sam- 
pling steps from the conditional distributions: 

„^CMB-Ei r CMB-S n n \ 

p\y h |c b ,9,c b ), 

p(c CMB-S |c CMB-E^ a)i 

} P (9\C™ B - E ,C™ B - s ,C b ). (29) 

The marginalized distribution for the CMB bandpow- 
ers, p(C & CMB - E ,C 6 CMB ~ s |Cb), with its associated covari- 
ance matrix, is then computed from the samples. This 
could be extended to include the SPT data, or data from 
Planck, for example. 

5.1.2. Calibration factors 

There are four ACT calibration factors. To minimize 
bin-to-bin correlations in the estimated CMB bandpow- 
ers due to calibration uncertainty, we divide out the 
148 GHz calibrations for the two ACT spectra, estimat- 
ing C' b CUB - E = C b CMB - E /yL for the ACT-E bandpow- 
ers, and C ' b CMB - s = C b CMB - s /2/L for ACT-S. 

We then estimate the 148 GHz calibration factors, 
j/i e , yis, and relative 218/148 GHz calibration factors, 
J/2e/yie, y2s/yis, as part of the secondary parameter set 
0. 

5.2. Marginalized CMB bandpowers 

Figure M shows the estimated CMB bandpowers from 
the ACT-E and ACT-S spectra, co-added together and 
compared to the multi-frequency spectra. The bandpow- 
ers for each region are reported in Table [5] Without 
assuming any cosmological model, the CMB bandpow- 
ers over the full angular range are remarkably consistent 
with the theoretical ACDM model predicted by WMAP. 
The uncertainty on the bandpowers rises at scales smaller 
than £ ~ 3000, and the correlations between bandpowers 
increases. 

Figure [TU] shows the effect of marginalization on the 
bandpower errors, using the ratio between the marginal- 
ized errors and the unmarginalized errors for a fixed sec- 
ondary model. It is clear that by measuring the spectrum 
at multiple frequencies, the CMB can be successfully 
separated from secondary contamination out to scales 
£ < 3500. At scales £ < 2000 there is little error infla- 
tion due to foreground uncertainty, and the errors are 
inflated by -20% (15%) by £ = 3200 for the ACT-E 
(ACT-S) spectra. The marginalized distributions for the 
CMB bandpowers are well approximated by Gaussians 
for multipoles to band-center £ = 3540 for ACT, as shown 
in the Appendix. 
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Fig. 10. — Inflation of errors due to foreground marginalization, 
relative to the errors for a best-fitting foreground model: at scales 
smaller than I ~ 2000, the errors are increased due to foreground 
uncertainty. 



We compare the secondary parameters recovered in 
this model-independent sampling to the case where 
ACDM is assumed. This comparison is shown in the 
Appendix; the parameters are consistent, with about 
a lcr shift in the estimated kSZ power. We find that 
the CMB bandpowers are not strongly correlated with 
the secondary parameters until scales well into the Silk 
damping tail at i }t 2500; a dominant correlation is then 
with the kSZ power due to its blackbody frequency de- 
pendence, and a smaller kSZ power - compensated by 
larger primary CMB power - is allowed when the ACDM 
assumption is relaxed. The CMB bandpower covariance 
matrix conserves this correlation information. 



TABLE 5 

LENSED CMB ANISOTROPY POWER 1 



l)C b /2-K (imK 2 ) 





ACT-E 


ACT-S 


Coadd b 


590 


2125 ± 156 


2286 ± 159 


2203 ± 111 


690 


1684 ± 113 


1662 ± 102 


1672 ± 75 


790 


2473 ± 143 


2209 ± 122 


2319 ±92 


890 


1922 ± 106 


1877 ± 100 


1898 ± 73 


990 


1007 ± 54 


1192 ±61 


1086 ± 40 


1090 


1214 ±62 


1124 ±57 


1165 ±42 


1190 


1028 ± 51 


990 ± 46 


1007 ± 34 


1290 


674 ± 33 


753 ± 37 


709 ± 25 


1390 


825 ± 40 


777 ± 35 


798 ± 26 


1490 


660 ± 31 


644 ± 29 


651 ±21 


1590 


445 ± 19 


476 ± 21 


458 ± 14 


1690 


387 ± 16 


357 ±17 


374 ± 11 


1790 


343 ± 14 


342 ± 15 


343 ± 10 


1890 


243 ± 11 


244 ± 12 


243 ± 8 


1990 


230 ± 10 


215 ± 11 


223 ± 7 


2090 


200 ±8.7 


206 ± 11 


202 ± 7 


2240 


138 ±4.9 


128 ±6.0 


135 ±4 


2440 


93.1 ±3.9 


88.9 ±5.2 


91.6 ±3.1 


2640 


58.0 ±3.3 


56.6 ±4.5 


57.5 ±2.6 


2840 


43.7 ±3.3 


45.6 ±4.3 


44.4 ±2.6 


3140 


23.2 ±2.3 


20.0 ±2.9 


22.0 ± 1.8 


3540 


9.8 ±2.4 


9.7 ± 2.9 


9.8 ± 1.9 



a To compute a likelihood using these data, the covariance ma- 
trix and bandpower window functions should be used, provided on 
LAMBDA. 

b This coadds the ACT-E and ACT-S estimated CMB bandpowers 
for plotting purposes only. 



5.3. The CMB-only likelihood 

We construct the CMB-only likelihood from the an- 
gular range where the CMB bandpowers are Gaussian, 
conservatively choosing I < 3500. We do not use the 
3500 < I < 4500 bandpowers as they are increasingly 
non-Gaussian, due to the foreground marginalization, 
and are more strongly correlated with foreground param- 
eters. The likelihood is given by 



21n^(C\ CMB |Cf) 



bidets. (30) 



Here 



CMB — E 



CMB-S 



Wbl.ACT-EjCf 1 
■ WM.ACT-sCf 1 



(31) 



where MB and X are the marginalized mean and co- 
variance matrix for the bandpowers, and Cf 1 is the lensed 
CMB spectrum generated from e.g., CAMB. We use 21 
bandpowers for ACT-E and ACT-S in the range 500 < 
I < 3500. A single calibration parameter fo r each region 
is mar ginalized over analytically, following Brid le et al.l 
(2002). The prescription for using this CMB-only like- 
lihood is simple, as no extra nuisance parameters are 
needed. 

To test the performance of this compressed likelihood, 
results are compared using both the full multi-frequency 
likelihood, and the CMB-only likelihood. Cosmological 
parameters are estimated for the restricted ACDM 6- 
parameter model, and a set of more extended models that 
probe the damping tail and peak shapes, including the 
running of the spectral index, the number of rclativistic 
degrees of freedom N e g, the lensing amplitude Al, and 
the variation in fine structure constant, a. Parameter 
constraints using both likelihoo ds agree to O.lcr, and are 
reported in lSievers et al.l (|2013| ) . We conclude that this is 
an efficient alternative to the full likelihood for the typ- 
ical extensions considered in cosmological analyses, al- 
though the full likelihood may give more optimal results 
for unusual models with features far into the damping 
tail. 

6. SUMMARY 

In this paper we have presented a likelihood formal- 
ism to describe the ACT multi-frequency power spectra 
that includes contributions from SZ and foreground com- 
ponents in addition to the lensed CMB. We model the 
data including four late-time astrophysical components: 
thermal and kinetic SZ, emission from CIB galaxies, and 
emission from radio galaxies. 

We have quantified these components using seven 
power spectra, splitting the CIB into a Poisson and 
clustered part, and including power from the cross- 
correlation between tSZ emission from clusters, and emis- 
sion from CIB galaxies that also trace the large scale 
structure. Rather than a minimal model with the fewest 
parameters, we have sought a model that includes our 
uncertainties with priors describing our knowledge from 
additional observations. For example, while the data do 
not demand that we include the tSZ-CIB correlation, we 
are motivated to include it to avoid placing unphysically 
strong limits on the kSZ power. 

Modeling these astrophysical components allows us to 
probe the primordial CMB fluctuations down to an angu- 
lar resolution of 4' using ACT. We have used the model 




Fig. 11. — Probability distributions of a selection of CMB bandpowers for ACT-E in the range 600 < £ < 4500, marginalized over 
secondary parameters (solid). The bin-center for each bandpower shown is indicated on each panel. The bandpowers are well-approximated 
by Gaussian distributions (dashed), except at scales I ^ 3700. The same behavior is seen for the ACT-S bandpowers. 



to extract an estimate of the primordial CMB spectrum 
well into the Silk damping tail, marginalizing over the 
foreground uncertainty. This produces a simplified com- 
pressed likelihood for use in cosmological parameter es- 
timation. 

We find that data observed by the South Pole Tele- 
scope give results consistent with ACT, accounting for 
the different removal of radio point sources, and differ- 
ent degree of contamination by Galactic cirrus. SPT 
and ACT have very different instrument design and scan 
strategies, and their observations on the sky have limited 
overlap. The excellent agreement between the datasets is 
not only an important cross-check but is another demon- 
stration of cosmic homogeneity. 
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APPENDIX 

In this Appendix, we perform additional tests on the CMB-only likelihood. In Fig [TT] we show a selection of the 
distributions of the CMB bandpowers from the estimated 600 < I < 4500 range. Distributions for the ACT-E and 
ACT-S bandpowers are compared to Gaussian distributions (dashed curves); bandpowers at £ > 3900 are significantly 
non-Gaussian, but are well fit by Gaussians at larger scales. The same behaviour is found for the ACT-S bandpowers. 

We then compare the secondary parameters estimated in two ways: (1) estimating CMB bandpowers, and (2) 
estimating 6 ACDM parameters. The distributions are shown in Fig[l2l and are consistent. The tSZ, point source 
parameters, and Galactic cirrus parameters are not strongly affected by the CMB model assumptions. The kSZ power, 
Oksz, is ~ lc lower in the model-independent case, as it is anti-correlated with the CMB bandpowers at I > 2000 due 
to the common blackbody dependence. The data cannot distinguish between lensed CMB power and kSZ power at 
£ ~ 3000 scales, so the preference for a smaller kSZ value in the model-indepedent case is driven by the prior that the 
CMB power is positive. We test this by allowing the CMB bandpowers to take unphysical negative values. Here, the 
kSZ power increases to aksz < 12 at 95% confidence, more consistent with the limits when ACDM is assumed. 
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